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A  class  of  data  windows  is  proposed  for  vtin  the  estimation  of  flu?  power 
spectra  cf  stationary  stochastic  processes.  These  windows,  which  are  genera  11  ca¬ 
tions  of  the  standard  Farzeu  filter,  me  constructed  by  the*  vise  of  appro,u  i ate ly 
norma li  zed  B-spl  i  nos .  It  is  demonstrated  how  this  window  c2r».;s  nay  be  computa¬ 
tionally  inipl  **,.»ent  ed  using  e  Fast  Fourier  Trans  torn  alnnri  thr.i,  inn  efficiency  cf 
the*  resulting  ptocedu.ro  is  general  Iv  a  significant  improvement  over  the  stale  of 
the  art  v<th  little  additional  computer  tine  required. 


INTRCWCf  fON 


TM  5  psper  fs  concerned  with  the  problem  of  the  result  of  usiii^,  lln.se  dr  la  windows  is  a  riv,- 

digitally  estimating  the  power  spectrum  of  a  vr$de-  nif leant  improvement  in  estimates  of  power  spectia 
sense  stationary,  ervodic,  stochastic  process  without  a  noticeable  increase  in  the  commit  At  i  c>nal 

from  a  sample  function  of  finite  length.  Not  effort.  Wo  have  chosen  to  me  continuous  notation 

long  after  Cooley  and  Tuhey  [5]  introduced  the  throughout;  conversion  to  discrete  formulae  where 

fast  Fourier  transform  algorithm  (FKTl ,  its  appli-  appropriate  is  straight  forward, 
cation  to  the  estimation  of  power  spectra  was  dls-  2  PREI IMINARIFS 

cussed  by  Bingham,  Godfrey,  and  Tukey  [2j  and  by 

Welch  [181,  They  suggested  a  method  based  on  com-  We  first  outline  some  basic  notation  and  re- 

plex  demodulation  and  pointed  out  the  computational  suits  which  are  assumed.  Here  we  follow  Stein  and 


tuts e  d r  l  a  v i ndews  i s  a  sig¬ 


nificant  improvement  ii'  estimates  of  power  spectis 
without  a  notlceehle  increase  in  the  commit  at  i  oral 
effort.  Wo  have  chosen  to  u?c  continuous  notation 
throughout;  conversion  to  discrete  formulae  where 
appropriate  is  straight  forward. 

2.  PRELIMINARIES 


We  first  outline  some  basic  notation  and  re¬ 


speed  advantages  which  can  be  obtained  by  using 
the  FFT  for  computing  Fourier  periodograms .  Welch 
proposed  a  direct  method  which  has  become  the 
standard  numerical  appioach  to  the  problem  of 
power  spectral  estimation.  This  procedure  has 
three  parts:  subdividing  the  available  data,  using 
^  ^  a  data  window  to  compute  a  smoothed  spectral  osti- 
9 ^  mate  for  each  segment,  and  averaging  these  smoothed 
|  spectral  estimates.  The  present  paper  proposes 
i  a  class  of  dafa  vitulww  derived  from  the  B-spllne 

basis  functions  of  Schoenberg  (14}.  f-.rrpi rlca I lv, 

fi,  . 


Weiss  (16}.  I,1  denotes  the  Banach  space  of 

Lebesgue  measurable  functions  which  arc  (lobes guc) 
integrable  over  the  real  line  F  .  For  fCb1  ,  the 
Fourier  transform  of  f  is  the  function  defined 
by 

hr)  =  r„  f(t)c'2"ictdt 

for  all  f,(F  where  i  =■  /-I  ,  The  notation 
(...)A  denotes  the  Fourier  transform  of  (...). 

For  ,  the  convolution  h  -■  f*g  is  defined 

as  the  function  b  C  Ll  given  by 
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for  .ill  ?  6  F.  If  f ,  g  r-  L  ,  then  (f*g)  -  fg. 

Similarly,  (fg)^  -  f*g  provided  Che  functions  in¬ 
volved  are  such  that  all  operations  make  sense. 

For  a  >0,  represents  a  dilation  operator, 

i.e,,  \f(t)  *  f(at).  Whenever  f,gc  L1 ,  it 

is  clear  that 

<‘.f*s.R>(n  =  l  M(f*g)C?)l 

and  0.f)*(?)  =  i*v  f(')  . 

We  shall  consider  a  sample  function  x(t) 

T  T 

on  a  finite  interval,  -  j  <  t  <  j  unless 
otherwise  stated,  of  a  stochastic  process 
( X(  t )  :  <  t  <  m) .  We  assume  this  underlying 

stochastic  process  to  be  wide-3ense  stationary 
(Doob{8|)  and  ergodic.  For  simplicity  we  take 
the  mean  value  to  be  zero,  bet  X(c)  have  the 
spectral  density  function  P(f).  The  terminology 
related  to  spectral  analysis  used  in  this  paper  is 
generally  that  of  Blackman  and  Titkey  *[  3 )  . 

it  is  apparent  that  high  resolution  and 
stability  (snail  bias  and  variance)  are  desired 
qualities  o f  puwe r  spec t  ra l  os t ima t os .  However , 
a  compromise  between  the  two  must  be  made.  Cn 
ordt.-r  to  obtain  a  stable  spectral  estimate  using 
the  classical  indirect  method,  an  appropriate 
function  of  lag  must  modify  the  sample  autcco- 
variancc  function.  Fundamentally  the  indirect 
method  estimates  P(f)  via 

e.  (f)  -  r  C(TlT.(  r)e'2rlfT  dr 

t  .  i:i  .  w  „  t 

whore  l.(T)  is  the  finitely  supported  lag  window 
and  cC"1)  is  the  sample  aulocovar i ance  function 
,  (T- It |)/2  . 

<:(T)  F  — (t-It|>/2  x(t+ 

Usually  the  support  of  I.(r)  is  much  less  than 
the  entire  interval  -T<  r<.  T  . 

The  at  .-aments  for  using  lag  windows  (e.g,, 
see  Jenkins  and  Watts  f  10)3  ran  be  ad-iptcd  to 
sliow  that  a  p:octdure  which  utilizes  the  FFT  to 
estimate  power  spectra  directly  should  incorporate 
a  data  window.  That  is,  for  a  direct  procedure 
to  have  a  proper  balance  between  resolution  and 
stability,  the  data  itself  must  be  modified  by  a 
;u l table  function  of  time.  Thus,  basically  the 


at  discrete  frequencies  where  W(t)  is  the 
finitely  supported  data  window. 

The  end  result  of  either  of  these  modifica¬ 
tions  is  that  we  in  fact  estimate  smoothed  values 
of  P(f).  More  specifically,  if  we  let  Pe(f) 
represent  the  estimate  computed  by  either  the  In¬ 
direct  or  the  direct  approach,  we  have 

*wg[P#(f)]=*  Q(f )*P(f ) 

where  the  average  may  be  taken  either  over  the 
ensemble  or  along  time  [3J.  Q(f)  is  called  the 
spectral  window  corresponding  to  the  window  used 
in  the  computation. 

The  direct  and  indirect  approaches  to  esti¬ 
mating  the  power  spectrum  of  a  stochastic  process 
are  not  funda.-.c'itally  different.  Blackman  and 
Tukev  [31  have  derive!  the  ccv.cc t i c*  between  a 
data  window  Wft)  and  its  t i valent"  lag  win¬ 
dow  L(~).  this  relationship  turns  out  to  he  the 
correlation  integral 

!.(->  -  ~  WfOW(fT)dt  . 

By  "equivalent"  we  mean  that  the  appropriate  use 
of  the  respective  "equiva l ant "  windows  yields 
power  spectral  estimates  which  are  equal  in  ex¬ 
pectation.  Note  that  if  the  data  window  is  an 
even  function  of  time,  then  the  equivalent  lag 
window  is  proportional  to  the  convolution  of  the 
data  window  with  itseLf. 

At  this  point  it  should  be  clear  that  in  the 
indirect  method,  the  spectral  window  Q(f)  cor¬ 
responding  to  the  lag  window  !.("■)  is  simply  the 
Fourier  transform  of  bfrl.  For  the  direct  method 
it  then  follows  that  t?»c  spectra  I  window  Q(f) 
corresponding  to  the  data  window  W(r  )  is  give*' 
bv  Q(f)  =  “lift))''  where  Iff)  is  defined  to 
he  the  Fourier  transform  ol  W(t).  J(fl  is 
called  the  frequency  window  correspondin'  to 
W(r).  We  again  refer  the  reader  to  Blackman  and 

Tukev  (31. 
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class  so  that  lh*' 


3,  srnr.T  data  v.-'hwws 
i.l  motivation 

lr  light  of  the  rat  ho  i  simple  relationship 
between  a  data  window  and  its  equivalent  lav  win¬ 
dow,  one  should  i;t>l  ignore  the  extensive  cifuit 
which  lias  boon  {ml  into  the  construction  of  lap 
window  hv  numerous  nut  hots  nor  the  proport i<s  of 
Lhose  windows.  Set*  for  exar.pJe  the  corp:  Jr.?  ir:i 
given  by  Jenkins  f  91  or  by  I'/trzen  [12,13). 

For  instance,  recall  the  continuity  class  to 
which  these  functions  belong.  In  parti  culm  we* 
observe  that  one  of  the  lag  windows  which  has  been 
proposed  by  Parzen  [13]  is  proportional  to  a  di - 
latod,  fourth  order  B-splinc  bRsis  function.  One 
of  the  two  data  windows  suggested  bv  Welch  (18] 
has  the. shape  of  a  second  order r  By-spline  basis 
function;  ho  recognized  that  his  proposed  window 
yields  tlu»  Parr.cn  spectral  window.  Henceforth 
we  will  refer  to  this  window  which  is  given  by 


as  the  Parzcn  data  window,  even  though  it  is  also 
a  member  of  the  family  of  spline  windows  proposed 
in  tin’s  paper.  Today,  Welch’s  procedure  (de¬ 
scribed  in  the  introduction)  incorporal i ng  the 
Par7.cn  window  is  undoubtedly  the  most  widely  used 
method  of  direct  power  spectral  estimation. 
Bingham  et  al  [2]  proposed  the  following  data 
window  which  is  continuous  through  the  first  de¬ 
rive  t  i  ve: 


It  I  <  AT 


W(t) 


.AT  <  It  I  <  AT 


N  0  ,  otherwise  . 

We  have  normalized  this  window  properly  although 
this  was  not  done  in  [2], 

From  another  point  of  view,  we  note  that 
under  the  Fourier  transform  operator  the  effect 
of  an  occurrence  at  one  point  in  time  tends  to  be 
spread  over  all  frequencies.  Furthermore.  time 
functions  with  corners  have  Four  ter  transforms 
that  have  a  "ringing11  effect.  Thus  It  seems  de¬ 
sirable  to  have  a  data  window  of  a  high  continuity 


function 

ore  to  e>:ar  :  v 
i  vs  bascM  upon 

A  pr el  in- 


class  ‘>o  that  lh*1  nvilifift  sampl* 
smoothly  to  /cm  at  it;:  nuls. 

Such  cons  i  not  nt  1  o-v^  r-.otivut. 
the  p  '^ibi  1  i  l  v  of  us**-’  data  \ni-i 
higher  ordei  r.-spliue  bps  it  !»?:.<  t  i 
iiiiii)  survey  is  ev.ro-.:  i  «v  i  ng.  Fiist,  we  i  ay  choc 
the  dr  si  red  continuity  class.  Sicondlv,  tin  coi  : 
poi.di::;*  spevt r*l  windows  have  th«  shape  ^  Hi — J 

so  that  the  magnitude  of  the  side  lobes  rapidlv 
decreases  relative  to  the  rain  pi  nr:  av  the  contin¬ 
uity  (or  order  k)  increases.  It  addition,  all 
spectral  windows  art:  of  r.cn-neget  i vt  type  (Parrel* 
f  1 5 1  ) .  l’hirdly,  dc'hoot  [7]  has  thrived  a  general 
algorithm  which  can  he  used  for  Accurately  evalu¬ 
ating  any  order  IJ-spline  basis  function  with  nnv 
i  lit  erval  %  of  support  and  at.  any  nurd*  length. 

3.2  THE  ITNoAMENiAl.  SliTN*  ivsn  loss 

Schoenberg  [ 1A ]  introduced  the  spline  funr- 
tions  in  19AA  with  the  following  definition.  V’ith 
k  a  positive  integer,  a  real  function  A(l)  de¬ 
fined  for  all  i  t  F  is  called  a  spline  f until  in 
of  order  k  or  degree  k-1  if  it  has  th.'  fol  If- - 
ing  properties: 

(O 

(U) 


(Ui) 


A(0  is  of  class  C* 

<A(t)  is  composed  of  polynomial  arcs  of 
degree  at  most  k-1  . 

The  polynomial  arcs  arc  joined  nowhere  other 
than  possibly  at  integers  n  If  k  is  even 
or  at  points  n4^  if  k  is  odd.  The 
spline  is  said  to  have  knots  at  these  points 
I f  we  define 


Jy-  . 

u 


y>  0 


l-i 


y  <  0 

with  exception  of  the  special  case  0+  ~  £  and 

let  6k  represent  the  order  central  differ¬ 

ence  operator  with  unity  step,  i.e.,  ft^fiy)  * 

f(y)  and  6kf(y)  -  bk~'  f(y+*  )  -  6k“1f(y-^) 
for  k  -  1,2,*..,  then  the  function 

M,(t>  -  (j^yr  eti'1 

is  easily  shown  to  he  a  simple  c:  ampin  oi  a  spline 

function  of  order  k,  M*(t)  is  called  the  funda- 


mental  B-splin.e  basis  function  of  order  k  .  Wo 
list  a  few  of  the  many  important  properties  of 
M*(t)  and  refer  the  reader  to  Schoenberg  [l'*,  15], 
Curry  and  Schoenberg  [6],  and  to  doBoor  [7j  for 

these  nr*J  other  properties.  M  v  ( t )  is  positive* 

k  k 

in  the  open  interval  (-  irt?r  )  and*  i dent i call y 
xero  elsewhere;  it  has  knots  at  -  j+v  where  v 
is  an  integer  satisfying  ^  <  v  <  k.  Mw(t)  is 
an  even  function,  and  its  integral  over  P  Is 
unity.  M*(t)  also  has  the  representation 


Mk(t)  * 


onrnr  £<-»> 


-!+, 


r- 


For  example,  the  fourth  order  function  is  given 
by 

r  i  /  .  or  .  ^  .  . 

iti  <  i 


K*(t)  = 


<  ■ 


i(2.|tl)>  -f(,- It  l)> 

iH't 


,  l<  It  I  <2 


.0  ,  otherwise. 

We  shall  require  two  integral  propertied: 


M«  (■)  '  j 


and 


CVM^Kt)  -  MpM(T)  . 

3.3  V\\K  PROPOSED  CLASS  OF  DATA  WINDOWS 
3.3.1  Ini roduct ion. 

The  class  of  spline  data  windows  we  propose 
consists  of  the  fundamental  spline  functions  each 
dilated  so  that  its  support  is  equa L  to  the  length 
of  the  sample  function  and  normalized  so  that  the 
total  power  of  die  data  remains  invariant  when  the 
spline  data  windows  are  applied.  If  we  let  Ck 
i  cpresent  tlu*  normalization  factor,  then  we  can 
use  the  results  of  .tin*  previous  sections  to  ob¬ 
tain  the  following  table  of  du*  spline  tint  a  w  in¬ 
flows  and  at*  1 1 » . :  corresponding  frequency  windows, 
lag  windows,  and  spectral  windows. 

Since  the  total  power  >r  is  given  by 

nr/ 

-  *  x*  (t)dt  *  T  Var[x|, 

*T/? 

we  seek  constants  C depending  upon  the  order 
k  of  the  spline  data  window  such  that  also 

a  -■  (cAeMOr'dt 
•  -V  3 


where  we  use  tlu*  lild.*  to  indicate  an  un -normal¬ 
ized  window.  The  expi-c  t  -*l  i  un  ot  this  equation  is 

<*  =  C-  '  ElV.dUKlx^d)  !<u 
'  2 


-  C;  Va  r ( k 1  '  W; ( t ) d  t . 

V‘T4 

Changing  the  variable  of  integration  we  obtain 
/  ^ 

c;  =  k/f  M?(t)dt  . 

i '  -+1 
■3 

An  application  of  Parseval's  theorem  yields 


0?  -  rk/l 


•ATh^r 


du . 


Therefore,  we  can  calculate  the  appropriate  nor¬ 
malizing  factor  for  any  order  spline  data  window: 
Cx  «  1.0  and 


C* 


2  VL-iliOL-ii-itlll 

.  i : t2k-i-i) : 


for  k''-^ . 


Normalizing  for  invariance  of  total  power  is 
equivalent  to  the  requirement 

1  (0)  -  j’"5  Q(f)df  -  l  . 


Data  Window 
Wk(t) 

C‘V«<U  '  (VK (a) 

Frequency 

W i ndow 

J»<0 

CVT  f i i n( -fT /k) V 

k  ’  \  ~ff  7k  j 

Equ  i  v.»I  ent 

Lag  Window 
U(r) 

T^VJ‘(r)] 

cS;  fk-r\ 

Spectral 

W  i  nd  uw 
%(l) 

i 

<1  /  J  in(  "fi/k)  \  2,1 

V.”  \  t{  l/~J 

TAB  IF.  1. 
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We  point  out  th;it  Wx  is  simply  a  trunca¬ 
tion  function,  l.j  is  the  lag  window  suggested 
hy  Bartlett  fl],  ^nd  is  the  correspond Lug 

liar t  Lett  spectral  window.  W3,  I,a ,  and  Q.,  are 
respectively  the  Farzon  data,  Lay,,  and  spectral 
windows . 

For  the  problem  of  balancing  i  lie  resolution 
and  stability  of  power  spectral  estimates,  ex¬ 
plicit  expressions  for  the  bias  and  variance  of 
smoothed  estimators  have  been  developed  in  the 
literature.  Jenkins  and  Watts  [10]  and  others 
give  these  expressions  for  the  indirect  computing 
method.  Welch  [17, 18]  derives  the  expected  value 
and  variance  of  the  discrete,  directly  computed 
estimates  for  a  stochastic  process  which  is 
tlauss i an ,  in  addition  to  our  previous  assumptions 
about  the  process.  In  each  case  some  property  of 
i he  window  being  used  is*  expl ie i t ly  involved, 
lb  a  practitioner,  however,  needs  seme  device  to 
aid  in  the  se l ec t  » on  of  a  v i e.dov  that  will  y i el d 
i  compromise  suitable  for  his  purposes,  Thu  most 
Cv-inon  measure  of  windows  is  c  tiled  bandwidth  (by 
inaJogy  with  filtering  p  rob  lens)  I 

Several  definitions  are  used  tor  bandwidth. 
Hackman  and  Tukey  [1]  use  inst.ad  the  term 
equivalent  width,  which  they  define  to  be  the 
ratio  of  the  square  of  the  Lute  •ml  of  tin*  spec¬ 
tra  t  window  to  the  integral  of  its  square.  if 
we  define 

l  P  (la(Odr, 

V  -30 

tin  n  the  equivalent  width  is  given  by 

»,  ('))  [[;"%(!  )>'• 

tonkins  anti  Watts  lloj  define  h  uidvidth  to  he 
th.e  width  of  the  rectnngul  «r  (••:»•, ertral)  window 
for  which  the  variance  ot  the  indirectlv  computed 
ni-»i  tril  estimate  is  the  same  .»■;  it  is  f  >r  t!ie 
a  i  ven  window.  I  hi  a  definition  of  bandwidth  vie  Id- 

I  I'M 

«,'<<?)  \  ■ 

A;t  Ins  boon  noted,  tor  a  properly  normalized 
window  WtOdt  1;  therefore  Ilf  .  We 


mention  BJ  for  its  application  to  f tie  spline 
lay  wind. «ws  given  in  lable  1.  However,  its  de¬ 
fining  ptopml/  is  no  long*  l  tun*  for  Lie*  direct 
method  of  computation  which  we  advocate.  For 
the  class  of  spline  windows  we  have 


Q?(fU’f 


•n:kvr  :  *-1  f-nWh'k  - 
k"  i !  1 4k-l-i)  * 

thus  Bj  (Q )  is  easily  determined  by  inverting 
ik  .  by  using  leas t -squares  teclmiques  for 
k  -  3 , , . . , 34 ,  we  obtain  the  approximation 

MQ*)  «.  (I.b8(k)°  *■*'*-  O.V2)/T  . 


Parzen  (t3|  defines  bandwidth  as  the  width  of  a 
rectangle  which  has  the  same  area  and  sane  maxi¬ 
mum  height  as  the  sp»«*tr»l  \vir.d*«.\  Ib.urc- 

fori-  '  *  ot  odr 

B.  *  0  »  - — 

max  Q(f) 
f 

and  l o r  the  spline  windows  we  have 


Another  coin -only  Us*l  :;VnsnV'’  of  -h'-slvi  dt  h  is*  the 
distance  *v  t  v#*«.  n  the  half-power  points  ot  the 
spectral  window,  lh.«;,  for  an  even  aptciml  win* 
dow  centered  at  t  0,  H,(Q)  ?t  where  T  is 

the  f  i  m  t:  frequency  such  that  0  (?)  1  Q  til).  r.et 

)n[*j  ’*•  ,  ■res-  :\t  1 1  •  *  sober  i  c  a  l  Hus., el  I  unctions 
of  the  tirst  kind,  then 
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variance  is  invetsel  pi  oport  i  onal  to  i  he  hand- 


.  O.th.  1‘  can  also  be  argu  d  that  t hr  bandwidth 
t  ;  1 1 •  window  riisi  bo  at  lea-it  as  small  as  tin* 
wi*'!!.  ui  the  Hat  i  uwv -!it  import- nl  detail  ia  the 

*  t  i  ii-i. 

iable  2  is  a  compilation  of  the  vaiious  con - 
s  la  its  associated  with  the  spline  windows  given 
ut  1  with  order  up  i  n  k-21  . 
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3.4  COMFIT  I  NO  WITH  SPLINT.  DATA  WINDOWS 
Except  fot  the  choice  of  data  windows,  the 
method  we  prop<>.;e  for  computing  power  spectra  is 
that  of  Welch  {18j.  Departing  from  our  conven¬ 
tion,  wo  assume  that  we  now  have  a  sample  <i*nc- 
lion  x(l)  ov ->r  the  interval  0<l^S.  Using 
for  example  a  bandwidth  criterion,  choose  the 
order  k  of  the  spline  window  to  be  used  and 
then  ^h<-  vfml'«»  length  7.  Subdivide  y(t)  inio 


N  segments  each  of  length  T  and  with  a  dis- 
pl.v.inciU  D  between  the  initial  points  of  con- 
setuiive  sepn«nts.  Denot  e  the  StgMnts  by  >;  ,  ( t  )  . 
Thus  (N  -1)1;  t  V  -  S,  and 

x;  (t )  -  x  1  (  j  -  I  >D  *t  ~  j  for  j  1 , . .  .  ,N;  f :  !  . 

Welch  uses  the  choice  of  D  unlv  .is  a  parameter 
to  increase  the  '.ability  of  the  resulting  isti- 
m-.ted  power  spec  Li  .aa,  In  the  c.it »  ol  splin<  win- 
do.;::  ;m  additional  consideration  should  be  tr.nae. 

I  1  D  -  f,  s oiue  o i  the  data  is  wo i  \ .h t  ed  substan¬ 
tially  heavier  than  the  remaindei .  As  can  be 
Seen  from  Figure  ],  the  inequity  is  more  pro¬ 
nounced  as  the  order  of  the  spline  window  in¬ 
creases  in  the  sense  that  the  spline  data  win¬ 
dows  become  progressively  more  peaked.  Since 
one  normally  assv.es  that  all  port  ions  of  his 
data  are  equally  representative  of  the  stochastic 
process,  intuitively  the  order  of  the  spline 
window  should  also  influence  I>,  In  practice  a 
good  rule  of  thivb  for  relatively  short  s.vsplc 

functions  seeirs  to  h°  h  =■  *T  for  k  -  4  I)  ^  t 

1  * 

foi  k  -  10,  and  D-  yT  for  k-  .'0  with  the 

ovei  lap  for  other  order  window/  t  a.*.».d  on  thi.  se. 

C 

Next  compute  the  smoothed  spectral  estimates 

H.(f)  «t|r^  x,(t)W.(t)<'*2T,if'dt  j*] . N. 

1  -  -t/S 

Finally,  the  estimated  power  spectrum  is  given  by 

r.(f)  -  s-  y-  Mf>  • 

w  J*  1 

Regardless  of  the  order  of  spline  window 
used,  P#(f)  is  obviously  a  consistent  estima¬ 
tor  of  P(f)  (at  least)  whenever  the  stochastic 
process  is  Gaussian.  'lo  see  this,  the  express¬ 
ions  for  Ell’*  (1)1  and  Var[l*#(f)]  as  given  by 
Welch  117,18]  need  only  be  compared  with  the 
convergence  properties  of  Liu*  spline  windows  ns 
discussed  in  this  paper.  Vnat  i ,  as  S  v*  ve 
have  the  following  type  convergence  for  a  sta¬ 
tionary  Gaussian  process:  if  T--* a  such  that 

1-0,  then  E(  (P(  f)  -  P, (O)a)-0  • 

In  implementing  this  procedure  on  a  digital 
computer,  one  could  derive  a  closed  form  for  the 

K-':jlino  basis  functions  like  that  gfvon  pie- 
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viously  for  M^ft).  However*  for  a  general  p^o- 
grnm  and  one  that,  would  undoubtedly  bo  numeri¬ 
cally  preferable,  the  algorithm  of  deftoor  (7)  is 
reciwwiiMidcd  for  evaluating  the  spline  functions. 

O f  c ou iso,  by  s yriino try  on i y  halt  th e  w i mi ow  ne etl 
be  directly  evaluated.  Brigham  (4}  thoroughly 
discusses  the  use  of  the  K'.'r  for  compnt  ing  the 
discrete  Fourier  transform  which  approximates 
the  Fourier  integral  we  used  in  writing  pj(£). 

4.  EXPERIMENTAL  RESULTS 

T‘he  effectiveness  of  the  spline  data  windows 
is  readily  apparent  from  experiments  with  machine 
generated  autoregressive  processes.  Power  spec¬ 
tral  estimates  computed  using  various  order  spline 
data  windows  are  compared  with  those  which  used 
the  B high am -God f rey-fukey  and  Parzen  data  windows 
and  with  estimates  computed  by  the  indirect  method 
jf  Blackman  and  Tukev  incorporating  the  Hanning 
lag  windi  j  [3).  The  estimates  were  nut  :;.al  ized  bv 
d  i vi  d lug  out  the  virl.ree  of  th.  tori-..;  Me  pro¬ 
cess;  as  is  customary  when  one-sidid  power  spectra 
a**e  used,  the  magnitude  has  been  doubled. 

For  these  experiments  we  used  sample  func¬ 
tions  of  -mi to r ogres s ive  processes  which  were 
machine  gene rated  according  to 


covariance  function. 

(ii)  B-G-*  :*v:  Hiuy, i.i  -(..nlf  rey-Tukey  «!.j»  a 

window  was  used  vi  th  segments  of  the 
sa..:)Je  function  of  length  r  v< 
samp l os . 

(iii)  Parven  xx;  Parzon  data  window  was  used 
with  segments  of  the  sat:. pie  function  of 
leii;;t!i  T  7  xx  s.tsip  I  os. 

(Lv)  Spline  xx-yv-zz:  Spline  data  window  of 
o i d o r  k  —  xx  was  us i*d  with  s e gr*.r : 1 1  s  of 
th«?  ;a::.ple  function  each  with  length 
T  -  yy  samples  and  displacement  D  -  zz 
samples  between  Miitiai  points. 

Mie  graphs  illustrating  these  experiments  arc 
found  in  Appendix  1. 


ojh  -t  1  ,  t  -*  l,...,luB4 


where  ::i  is  the  order  <>f  the  process,  (*vl }  are 
the  regression  coefficients,  and  e,  is  a  white 
noise  process  with  E|-r.  j  U  and  Var(s,  )  -  1. 


i  in’  parti*  ular  in  oc*. 


.vaich  were  used  in  the 


i  i  !  u  ;  l  r.F  ions  had  the  following  coot  f  j  c  ients : 

(i)  First  order,  -0.40 

(ii)  Second  order,  ..  1.00,  -0.40 

(iii)  Fourth  order,  ^  0.73,  0.40, 

,,  -  -0.33,  ...  *  -0.40  . 

Ihe  following  seller;.'  identifies  the  windows  in- 
vu!  /-'I  In  i  !m  c  .-input  at.  jam  of  the*  various  power 
spectral  e.ti mates, 

(i>  dinning  xx :  estimate  computed  by  the  in¬ 
direct  tFod  with  the  Hanoi nz  lag  window 
[1]  and  using  x<  lags  of  the  simple  auto- 
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